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Abstract 



7 ! ' Mills, Robbins, and Rumsey [g conjectured, and Zeilberger ]1^ recently proved, that there 

^j. . are ^u'^Wy} Aw-i)! alternating sign matrices of order n. We give a new proof of this result 

PsJ I using an analysis of the six- vertex state model (also called square ice) based on the Yang-Baxter 
equation. 

\m/ , Mills, Robbins, and Rumsey P conjectured that: 

(-i I Theorem 1 (Zeilberger) There are 

"t^ ■ 

a- ^^ ^ l!4!7!...(3n-2)! 
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n!(n-hl)!(n-F2)!...(2n-l)! 
n X n alternating sign matrices. 

^ ' 

1^^ ' Here, an alternating sign matrix or ^S'M is a matrix of O's, I's, and — I's such that the non- 

f^ . zero elements in each row and column alternate between 1 and —1 and begin and end with 1, for 

^N ' example: 

(N 

-(— > 

2 I Alternating sign matrices are related to a number of other combinatorial objects that, remarkably, 

H ■ are also enumerated or conjectured to be enumerated by ratios of progressions of factorials or 

^ . staggered factorials p, . 

k>( I Zeilberger [^] recently proved Theorem |^ by establishing that ASM's are equinumerous with 

Vh ' totally symmetric, self-complementary plane partitions, which were enumerated by Andrews [l[. 

In this paper, we present a new proof. The most interesting part of the proof is due to Izergin and 

Korepin |^, ^], who follow Baxter's remarkable use of the Yang-Baxter equation Q. 

If X is a number, define the x-enumeration A(n; x) of n x n ASM's as their total weight, where 

the weight of an individual matrix is x^ if it has k entries equal to —1. A variation of the proof 

establishes another conjecture of Mills, Robbins, and Rumsey: 

Theorem 2 ASM's are 3- enumerated by 

^ V {n + iy.{n + 2)\ . . . {2n)l 



A second variation establishes the weU-known 2-enunieration of ASM's [^, ^]: 

Finally, the following result, also conjectured by Mills, Robbins, and Rumsey, follows easily from 
the general method: 

Theorem 3 For each n, there exists a polynomial B(n;x) such that 

A{n; x) = B{n; x)B{n + 1; x) 

for n odd and 



A{n; x) = 2B{n; x)B{n + 1; x) 



for n even. 



Mills, Robbins, and Rumsey further conjectured that for n odd, B{n;x) is the x-enumeration 
of vertically symmetric ASM's (where the weight is x'^ if there are k ones to the left of the middle 
column), but this relation remains open. 
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1 State sums 

The six-vertex model in general refers to the multiplicative weighted enumeration of orientations 
of a tetravalent planar graph G (called states) such that at each vertex, two arrows go in and two 
go out. Number the six allowed orientations incident to a given vertex (called states of a vertex) 1 
through 6: 

1 



1 



State i at vertex v is given a weight w{i,v). The weight of a state of G is the product of the weights 
of its vertices, and the state sum is the total weight of all states. The six-vertex model may also be 
considered with boundary conditions, meaning that there may be univalent vertices whose edges 
have fixed orientations. In particular, consider a six- vertex state of an n x n square grid with edges 



pointing inward at the sides and outward at the top and bottom: 
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The six-vertex model on a square grid is also called square ice. A square ice state can be converted 
to an ASM by the correspondence 



















This conversion is bijective |^, ^]. Thus, the enumeration of ASM's is equivalent to a six- vertex 
state sum in which all weights are 1. 



Let hhe a complex number or an indeterminate, let q^ denote e , and let [x] denote ^ 



7^/2_ 



-x/2 



q'-T^-q- 



1/2 



We will consider various half-integral Laurent polynomials, meaning polynomials with integral or 
half-integral exponents of either sign such that the difference between any two exponents is an 
integer. For example, if q is fixed, [x] is a half-integral Laurent polynomial in q^ . Given two such 
polynomials P{t) and Q{t) over a ring A, we will say that Q divides P if P{t)/Q(t) G ^[t-^/^, t~^/^]. 
For example, t divides 1. 
A vertex labelled by x: 




-x/2 



x/2 



[x-l] 



[x-l\ 



(1) 



\x\ 



(Since the weights are invariant under rotation by 180 degres, but not 90 degrees, the meaning of a 
vertex depends on which pair of kitty-corner quadrants contains its label.) Such a vertex is called 
an i?-matrix and is also denoted as R{x). 



Theorem 4 (Baxter) If x = y + z, the R-matrices R{x), R{y), and R{z) satisfy the equation 





This remarkable identity is known as the star-triangle relation or the Yang-Baxter equation|2|. 
Specifically, R{x) is said to parameterize the trigonometric solutions to the Yang-Baxter equation. 
Before proving it, we discuss exactly what the equation means. Each of the two graphs in the 
equation has six external edges, meaning edges with a univalent vertex. For each external edge on 
the left, there is a corresponding external edge on the right whose univalent vertex is in the same 
position; for example, on both sides there is a lowest univalent endpoint, and the two edges with 
this endpoint correspond to each other. For each of the 64 orientations of the external edges on 
the left, one can form a state sum Z by summing over admissible orientations of the three internal 
edges, and one can consider the same orientation on the right and form another state sum Z' . The 
equation then says that the Z = Z' m all 64 cases. In order for the state sum to be non-zero, three 
edges must point in and three must point out, so the identity is trivial in 44 of the 64 cases. Note 
further that the equation simply says that the left side is invariant under rotation by 180 degrees, 
so the 20 non-trivial numerical identities reduce to 10 identities repeated twice. The argument 
that follows uses other tricks to further reduce the number of numerical identities to one which can 
checked be checked easily: 

Proof: We first rearrange the left side of the Yang-Baxter equation: 




Consider the following augmentation of the six-vertex model: Suppose that a graph has a curved 
edge with a horizontal tangent at a point p and which is concave down at p. If the edge is oriented 
to the left in some six- vertex state, p is assigned a multiplicative weight of —q^''^, but if it points to 
the right, it is assigned a multiplicative weight of 1. Contrariwise, if the tangent is horizontal but 
the curve is concave up, p has weight — g~^'^ when the edge points to the left and weight 1 when 
it points to the right. With this convention, the following simple identities hold: 




.ql/2 _ g-l/2 



-[2] (2) 



Moreover, R{x) can be expressed as 





+ [x-l] 



Thus, a six- vertex state sum involving i?-matrices can be expanded as a sum of curves in a calculus 
in which each closed loop contributes a factor of —[2]. (This calculus is called the Temperley-Lieb 
category and is closely related to the quantum group C/g(sl(2)) [^.) The calculus is invariant under 
isotopy of curves by equation (|2[). The left side of the Yang-Baxter equation then expands to 
eight terms, which may collected into five terms corresponding to the five crossingless matchings 
of six points on a circle. Three of the matchings are invariant under rotation by 180 degrees. The 
coefficients of the other two are 

[z-l][x][y-l] 

and 

[z][x -l][y-l] + [z][x][y] + [z-l][x- l][y] - [2][z][x - l][y]. 

These two quantities are rendered equal by the identities x = y + z, [—a] = —[a], and 

[a][b] - [a + l][b - 1] = [a - b + 1]. 



Thus, the left side is invariant under rotation by 180 degrees. 
As a final notational convenience, define 



D 



x-y 



y 



when the lines rather than the vertices of a tetravalent graph are labelled. Following Izergin 
and Korepin |5|, ^, consider nx n square ice with arbitrary parameters X = xo,...,a;„_i and 
Y = yo, . . . , y-n-i for the horizontal and vertical lines: 
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Xn-1 — >- 



yo yi 
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Vn-l 



Let Z(n; X, Y) be the resulting state sum. 

Lemma 5 (Baxter) The function Z{n; X, Y) is symmetric in the Xi 's and in the yi 's. 

Proof: Consider the iih. and i + 1st horizontal lines. An extra vertex (implicitly labelled by 
Xi — rcj+i) may be introduced on the left at the expense of a generically non-zero multiplicative 
factor: 



Xi+l 



Xi+l 



[Xi - Xi+l - 1] 





This relation holds because in an allowed state, all four edges of the new vertex must point to 
the right. By the Yang-Baxter equation, the vertex can be moved from the left side to the right, 
whereupon it can be removed, which recovers the multiplicative factor. This operation switches 
the labels Xj and Xj+i. Therefore Z{n;X,Y) is symmetric in Xi and Xj+i for each i, which renders 
it symmetric in all Xj's. The same argument applies to the yj's. □ 

Lemma 6 If Xi = yj + I, then 

Z{n;X,Y) = -q-'/^\\\[xi-yk] I I \\[xk-yj]\Z{n-l;X\xi,Y\yj). 



Proof: Assume first that i = j = 0. By Figure (^), the upper left vertex must have state 1 in 
a non-zero state of the grid. This forces the rest of the top row to have state 5 and the rest of 
the left column to have state 6, which yields the given multiplicative factor. (In terms of ASM's, 
only those matrices with a 1 in the top left corner contribute.) The remainder of the grid is an 
n — 1 X n — 1 square ice state. 

The general case follows from Lemma ^. □ 

Lemma 7 The quantity q"'^'-''''^Z{n;X,Y) is a polynomial in q^*^ of degree at most n — 1. 

Proof: If we multiply all weights of vertices in the first row by 5^0/2^ then q^° appears linearly 
in those weights in which it appears at all. Therefore the modified state sum 

Z'(n) = g"^°/2^(n;X,y) 

is a polynomial in q^". The first row is the only row in which xq appears. In this row, there must 
be one vertex in state 1, whose modified weight does not involve xq, and n — 1 vertices in state 5 
or 3. (In terms of ASM's, there must be a 1 in the top row.) Therefore Z'(n) has degree at most 
n- 1. □ 

Theorem 8 (Izergin,Korepin) The state sum Z{n; X, Y) is given by 



(-1)'" nr=o' q^y^-^^^/') Uo<.,<n[^^ - yj][^^ - y, - 1] 

Z{n; X, Y) = j) \7 \ ^^* ^' 

[Y\o<J<^<n[^i - Xj]) [Uo<i<j<n[yi " ^il^ 



where 



''^ [xi - yj] [xi - yj - 1] ■ 



Proof: Lemmas and ^, together with ^(0) = 1, inductively determine Z{n) by Lagrange 
interpolation. It is routine to check that the right side satisfies Lemma ^. To check that it also 
satisfies Lemma ^, Let P be the numerator, let Q be the denominator, let D be the determinant, 
and let D' be a term in the expansion of the determinant. The product PD is a half-integral 
Laurent polynomial because PD' is for any choice of D' . Moreover, Q divides PD, because D is 
antisymmetric in the Xj's and in the y^'s and therefore in the g^*'s and in the g^^'s. Thus, qD 
is a half-integral Laurent polynomial polynomial in g^o. Finally, the leading term (expanded as 
a Laurent polynomial in g^") of any PD' has exponent (2n — 3)/2, while the trailing term has 
exponent (1 — 2n)/2. Therefore the same is true of PD, and qD has leading exponent at most 
(n — 2)/2 and trailing exponent at least — n/2. In conclusion, q'^^^'^pD is a polynomial in g^" and 
has degree at most n — 1. □ 



2 Determinants 

Consider the state-sum value 

Zi(n) = Z(n;^,^,...,i,0,0,...,0). 

In any n x n square ice state, there are n more vertices in state 1 than state 2, equal numbers in 
states 3 and 4, and equal numbers in states 5 and 6. Since the weights of these states in R{2) are 
— g~^'^, —q^'^, — [|], —[^], [^], and [i], respectively, it follows that 

A(n;x) = [i]"-"'(-l)V/^Z.(n), (3) 

where x = l/[2]^ = [2] + 2. Unfortunately, the determinant in Theorem |8| is singular for Zi/2(ra)- 
Therefore, we will instead evaluate 

Zi (n; e) = Z{n; - + e, - + 2e, . . . , - + ne, 0, -e, -2e . . . , (1 - n)e) 

when h = ^g , which implies that x = 1. 

Let s = q^. Firstly, 

lu, li s'^ + l + s-'' 

Ike + - Ike = 

^ 2^^ 2^ -3 

and 

k/2 _ -k/2 

\ke] = ^— . 



The matrix M of Theorem H becomes 



M, 



The state sum becomes 

^-n/4^-„V23-n(n+l)/2 j-[^^^ .<„(s^+J'+^ + 1 + s-(^+^^+l)) 

^^ (^' ^) = n.. ■ (,sC^-i)/2 - .0---)/2)2 det M. 



The determinant of M can be computed using the following two lemmas, for which we extend the 

-l-X /2 7- X /2 

bracket notation by defining [x]t = — ~ i — for any t. 

Lemma 9 (Cauchy) Let X = xq, . . . , 2;„_i and Y = i/q, . . . , Hn-i be variables, and let 

1 



Tin;X,Y) 



''^ [xi - yj]t 
for < i,j < n. Then 

,i-i-0<j<i<niXi ~ Xj\t) {[[o<i<j<niyi — VjU 



det T(n, k;t) 



[[o<i.i<niXi yjU 



Proof: Let D be the determinant, let P be the denominator, and let Q be the numerator. Then 
the arguments of Theorem g apply, but with the conclusion that qD is a, degree polynomial in 
all variables, i.e., a constant. Let D' be the diagonal term in the determinant; D' is the only term 
such that PD' is not divisible by any [xi — yi\t. All factors of qD' cancel at the specialization 

Xi = Hi and all other terms of qD vanish; therefore qD = 1. □ 

Lemma ^ can also be proved by induction using Dodgson's condensation method y, IC]. 
Let T{n) = T{n; 1, 2, . . . , n, 0, -1, -2, . . . , 1 - n). Then 

T{n)i 1 = — 

and det T{n) is given by Lemma |9|. 
Lemma 10 Let 

g{i+j+l)/2 _ ^-(i+j+l)/2 

for < i,j < n. Then 

(_l)n(n-l)/2 / .^ \ / .^ gl/2^(i-i)/2 _ gl/2^(j-i)/2 

det 5'(n;s,t) - ^ ' I 1 r r. .i2 i i i r 




(tl/2_t_l/2)n 1 ii ' '' \ -'■-'- [i+j + l] 



Proof: The quantity s" '^(det 5(n; s, t) has degree re^ as a polynomial in s. Moreover, for 
<k <n, 

k—l 

.i-{k~l)/2^ 






S{n;t'',t)ij = [ky+j+i =J2A{t 

where the matrix A{z) given by 

A{z\, = z^+^+' 

has rank 1. Thus, the rank of S{n;t^,t) is at most k, and it follows that (s — t^^"-~^ divides 

det ^(n; s, t). Similarly, (s — t~ )^~ divides the determinant for 1 < k < n. These divisibilities 

determine detS'(n;s,f) up to a factor which is a function of t. The leading coefficient is then 
(detT(n))/(tV2_t-i/2)n_ □ 

The determinant of —M/3 = S{n; s, s^) is given by Lemma ^. Collecting factors yields 

Zi/2{n;e) 



^_„/4^-„V23-n(n+l)/2 no<.,,< J3(^ + j + 1)] 



no<^,,<n[i+j + Ms)(n 



0<j<i<n 



(_3)n(_i)n(n-l)/2 ^ TT i^i ' j)?s\ I TT [3].[3(i - j) + 1], 




j=o nj=ib]s / 

Note that the second factor is Andrews' g-enumeration of descending plane partitions |9[ with q 
replaced by s. Taking the limit as e — > and combining with equation (^), the factors of q and —1 
cancel, the factors of s become factors of 1, and the brackets disappear. The result is 

l!4!7!...(3n-2)! 
n!(n + l)!(n + 2)!...(2n-l)!. 



This completes the proof of Theorem |^. 
For general x, the matrix M becomes 



M, 



x^ — 4x 



There are two other values of x when the denominator is a cyclotomic (Laurent) polynomial in some 
power of s, namely x = 2 and x = 3. In the former case, — M/4 = S{n\ s^, s^), whose determinant 
is given by Lemma |lO|; alternatively, the determinant may also be derived from Lemma p. In the 
latter case, — M/3 = S"(n;s,s^), where 

S (n; s, t)ij = --igipj-pY. 

A variation of Lemma |l^ establishes the determinant of S'{n; s, t); the leading coefficient is simply 
the determinant of S{n;t,t^). These manipulations clearly lead to product formulas for A{n]2) 
and A{n;'i), and in particular, to a proof of Theorem g. We omit the details of rearranging and 
cancelling factors to put the product formulas in their standard form. 

Finally, we use Theorem ^ to prove Theorem y. Recall the variables Xi and yi in the definition 
of Z, which are not to be confused with the x of A{n; x). If we set Xj = ^ + fie and yj = fje for 
some /j's such that /„_i„j = — /j, then Z{n]X,Y) again converges to A{n]x) up to normalization 
as e ^ 0. In this case, the corresponding matrix M is given by 

x^ — Ax 



This matrix M possesses the symmetry (i, j) i-^ (n — 1 — z,n — 1 — j), i.e., it commutes with the 
antidiagonal permutation matrix P. Therefore, a change of basis divides M into blocks correspond- 
ing to the eigenspaces of P. Therefore the determinant of M is the product of the determinants of 
the blocks. This is the origin of the factorization of the A{n] x)'s into the B{n] x)'s. 
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